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Abstract 



Constructive techniques for controlling a coupled, heteronuclear spin system, via bounded 
amplitude sinusoidal pulses are presented. The technique prepares exactly any desired 
unitary generator in a rotating frame, through constant controls. Passage to the original 
coordinates provides a procedure to prepare arbitrary unitary generators, via bounded am- 
plitude piecewise sinusoidal pulses whose frequency is one of the two Larmor frequencies 
and whose phase takes one of two values. The techniques are based on a certain Cartan 
decomposition of S'f/(4) available in the literature. A method for determining the pa- 
rameters entering this Cartan decomposition, in terms of the entries of the target unitary 
generator, is also provided. 



1 Introduction 

The goal of this paper is to provide, constructively and exactly (i.e., without any approx- 
imations) a decomposition: 

satisfying i) Ol > 0; and ii) 02 \hk |< C*, for some prescribed bound C . The matrices, 
lij^i = l,2,j = x,y,z, defined in Section 2, are the standard tensor products of I2 and 
the Pauli matrices. As will become clear, later the same decomposition can be refined to 
simultaneously satisfy \ — \< D for another bound D. This problem arises in the control 
of a coupled, heteronuclear, two spin system which is being controlled by addressing each 
spin individually. The adjective "heteronuclear" is , throughout this paper, meant only to 
signify that the Larmor frequencies of the two spins are different. The condition Ol just 
means that the time for which a pulse has to be applied must necessarily positive, while 

02 (or more precisely, the condition | ^ |< D) represents the fact that hard pulses are 
not being used. 

Coupled spin systems are useful from several points of view. They are certainly ubiq- 
uitous in NMR studies, ^ ^. They also provide examples of coupled qubits in quantum 
computation and information processing, ^ ^ |6j . One method to control spin systems 
is to use hard pulses, i.e., high amplitude pulses applied for very short times. The usage 
of hard pulses, though useful from the perspective of minimizing the time consumed to 
prepare a desired unitary generator, 5, has its own problems. Perhaps most importantly, 
the usage of hard pulses allows one to only approximately prepare S, since the available 
resources in any situation is limited. For certain problems of NMR spectroscopy this may 
not be such a big restriction since the translation of a certain NMR objective to the prob- 
lem of unitary generator preparation may allow for a broad choice of unitary generators 
(though, to the best of our knowledge, there has been no systematic assessment of the 
inaccuracies introduced by the usage of hard pulses for specific NMR objectives). How- 
ever, for purported quantum computation applications it is indeed a restriction, since the 
desired S has to be obtained with great accuracy. Secondly, the usage of hard pulses may 
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violate the basic feature assumed often in the methodology - namely the ability to address 
single spins selectively. This already necessitates restriction to heteronuclear molecules. 
However, even for heteronuclear molecules the usage of infinite amplitude is problematic. 
Indeed, neglecting modes which are coupled to a system being studied is essentially a 
perturbation theory argument. Thus, not only has the frequency of the pulse got to be in 
resonance with the subsystem being studied, but the pulse area has also got to be bounded. 
If high amplitudes are being used then the pulse has to be applied for an extremely short 
time. It is unrealisitic to assume that a pulse can be applied for an infinitesimal time 
in the laboratory. The situation is analogous to the usage of the rotating wave approx- 
imation in molecular control studies, 0, ||, Not only must there be no resonances 
amongst the coupled levels, but the amplitude of the pulse must be much smaller than the 
frequency separation between the pair of levels being addressed. Similar considerations 
occur in NMR spectroscopy, pc| ]. These problems are compounded further in quantum 
computation applications, since it is desirable to perform local operations on individual 
(or small collections of) qubits, without too much crosstalk with other qubits during these 
operations. Usage of hard pulses, will eventually cause other qubits to be coupled. Put 
differently, one of the goals of hard pulse technology, namely the avoidance of decoherence 
effects, can in fact be defeated by the very usage of hard pulses. Thus, while it is desirable 
to finish all control action before relaxation processes become dominant, it is even more 
important to use fields which do not cause other couplings or processes neglected in a 
model to become significant. Quite often this means essentially that the field be bounded 
in amplitude. In this paper, arbitrary bounds on the amplitude will be allowed. Thus, while 
uji 7^ UJ2 (where the LUi,i — 1,2 are the individual frequencies) is needed, the difference need 
not be as large, as would be required by the usage of hard pulses. Thus, there is lots to 
be said for controlling spin systems without hard pulses. In this paper, explicit techniques 
which avoid hard pulses are provided for a coupled spin system. 

It is worthwhile to place in proper mathematical context the question being studied. 
Coupled two spin systems can be viewed as examples of left-invariant systems, with drift, 
evolving on the compact Lie group, SU{4:), ||l^. It is known non-constructively that, 
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under a certain Lie algebraic condition, such systems can be controlled with piecewise 
constant pulses whose amplitude can be arbitrarily bounded [ pl . Providing constructive 
proofs for such results, on the other hand, is an entirely different matter. There are 
two sources of complication. First, the ambient Lie group, S'C/(4), is high-dimensional 
(precisely, fifteen dimensional). But more importantly, the presence of a drift term (i.e., 
the free Hamiltonian) significantly complicates constructive control. Indeed, from a control 
theoretic perspective hard pulse arguments are ways of avoiding the effects of drift (though 
at the expense of introducing far more deleterious effects) . There is an extensive literature 
on constructive control for driftless systems - see, for instance, the survey, [l^ . The same 
reference contains a survey of special classes of classical mechanical systems with drift, 
which can be controlled constructively. In [|l3j, building on earlier work a detailed 
study of constructive control of systems with drift on SU{2) was provided. The special 
structure of SU{2) played a crucial role in this effort. In this paper, we will show that the 
special structure of SU(2) allows, once again, to constructively control the spin system 
being studied here. This structure of SU{2) enters in two manners. First, the groups 
SU{4) and 5*0(4) are intimately connected to one another due to one well-known Cartan 
decomposition of SU{4). This Cartan decomposition factors every S £ SU{4:) as products 
of matrices in SO (A) and matrices which are exponentials of a Cartan subalgebra of the 



complement of so(4), the Lie algebra of S'0(4), in su(4), the Lie algebra of S'[/(4) [|_4, 15 
[T^ . Now 80(4) is essentially the same as SU{2) (g) SU{2). Not only does this fact clearly 
establish this Cartan decomposition (this is well known), but it also enables techniques 
inspired by our earlier paper, for systems on SU{2) to achieve the main goals of this 
paper. The second role of SU{2) is to facilitate the calculation of the real parameters in 
this Cartan decomposition. In the appendix, a method to determine these parameters as 
explicitly as possible (explicit modulo the solution of a pair of transcendental equations) 
is provided, and once again the structure of SU(2) is the main ingredient. 

In the authors use this particular Cartan decomposition of SU{4) to address the 
problem of generating any unitary generator S G SU{4:) for a coupled two particle spin 
system via selective one-spin hard pulses. The essential difference between the approach 



taken here and that in |17[ | stems from a Lie theoretic nuance and may be summarized as 
follows. This Cartan decomposition and plus a few calculations leads to a factorization of 
any S G SU{i), 0: 

S = n'^^^expi-itkAk) (1.2) 
where each Ak is one of IizI2zt Iix, I2xt hyi l2y Greater detail about the decomposition, 



(1.2), is provided in Equation (2.10) in the next section. Under the assumption that it is 
possible to address single spins selectively, (and a subsequent passage to a rotating frame) 
Iizl2z represents the internal Hamiltonian, whereas the remaining matrices represent the 
control coupling. In factors of the type Qxp{—itt.Iizl2z) are generated by free evolution, 
whereas factors which are exponentials of the control couplings are generated by hard 
pulses, so that the drift term, i.e., Iizl2z, makes a negligible contribution to such a factor. 
This, naturally, makes the generation of the desired S (even in the rotating frame) at best 
approximate. In our approach, we will also generate the factors exp{~itklizl2z) through 
free evolution. However, instead of viewing the remaining factors as the exponentials of 
the control couplings, we view them as the exponentials of the iterated commutator of the 
internal Hamiltonian and the control couplings. In this iterated commutator, the internal 
Hamiltonian occurs twice whereas the control couplings occur only once. Thus it may 
be surmised, based on our prior experience p3[ , that such exponentials can be generated 
via three factors, two of which are free evolution terms and one is a control (though not 
hard) pulse. While, this conjecture has not been established for general nonlinear control 
systems, it turns out to be valid for the two spin system studied in this paper. Put 
differently, in generating the exponentials of the control couplings we make systematic 
use of the drift term, as opposed to treating it as a nuisance which can be overcome via 
hard pulses. As a further step towards making the whole procedure constructive methods 



which yield the tk in Equation (1.2 ) in terms of the entries of the target unitary generator, 
S, will be displayed. The fact that the Cartan decomposition is the analogue of the 
Euler decomposition for the higher dimensional unitary groups is well known, ||lj, |l5|, ^ . 
However, a stumbling block in its usage is the fact that the analogues of the Euler angles, 
i.e., the real parameters tk, have never been explicitly written down in terms of the entries 



of the target matrix S. 

Ultimately, the choice of hard versus soft pulses depends on both the system and the 
target in question. The contribution of this paper is to show that many objectives can 
be met via soft pulses. Furthermore, this is demonstrated by making use of the iterated 
commutators of the internal Hamiltonian and the control Hamiltonians, in contrast to 
hard pulse approaches which make no use of this commutator. This is another reason why 
the results below may be of interest. Indeed, all arguments for the controllability of finite 
dimensional quantum systems |l9|, ^ ^ hinge on the commutators of the 

internal Hamiltonian with the control Hamiltonians. Further, these arguments also show 
that it is possible (non-constructively) to prepare any target via bounded amplitude fields. 
Thus, it is at least didactically pleasing to demonstrate controllability constructively via 
methods which explicitly use such iterated commutators. In summary, two important 
problems are solved constructively in this paper: a) the control of a coupled spin system 



via sinusoidal, bounded amplitude pulses assuming the factors tk of Equation (1.2) are 
given, and b) finding these parameters, tk, algorithmically (so that there is something 
useful even for afficionados of hard pulses). For both these problems a significant role is 
played by the structure of SU{2). 

The balance of this paper is organized as follows. In the next section some notation and 
basic facts about the Lie group SU{2) will be collected. In particular, an explicit formula 
for an Euler angle factorization, with factors exponentials of ax and (jy , is provided. The 
third section carefully derives the rotating frame for the basic model, and shows how 
to determine the frequencies and phases of the piecewise sinusoidal pulses which will be 
designed to generate any S G SU{4). The fourth section describes how to determine the 
amplitudes of these pulses and in the process explains how the two goals, Ol and 02 
mentioned at the very begining of this section, help in finding the amplitude and the 
duration of the pulses. Some conclusions are offered in the next section. An appendix 
describes how to go about calculating the Euler angle analogues for S'C/(4). 
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2 Review of Basic SU{2) Facts 



The Pauli matrices will be denoted as: ax 



^ ^ 





i 



and (7» 



v 



^1 Oj 

In terms of this definition of the Pauli matrices, a few important 4x4 



-1 ^ 

Hermitian matrices can also be defined: 



hk ^ cTk® h]hk ^ h®cFk, k = x, y, z; 



(2.3) 



Note that both in the definition of the Pauli matrices and in the above equation, the 
customary factor of i has been omitted. This is for notational convenience and does not 
affect any of the results below. 

Next, a very useful representation of the SU (2) matrices, the Cayley-Klein parametriza- 
tion, follows. 



S = S'(a,C,M) = 



e"" cos a e'*^ sin a 



V 



(2.4) 



a, C, ^ are the Cayley-Klein parameters of S. Since this parametrization is nothing but 
the entries of S written in polar form, it is clear that ^ and /i may be taken to be in [0, 27r) 
and a to be in [0, §]. 

Crucial for the purposes of this paper is an Euler parametrization, in terms of and 
(7y, of any S'(a,C,/i): 



(2.5) 



Though, in principle, the Euler angles D, E and F may be obtained from more common 
Euler angle parametrizations p3| , an explicit formula for them in terms of the Cayley- 
Klein coordinates goes a long way towards making the techniques of this paper genuinely 



constructive. Therefore, the following relations are very useful, |24 



cos(£;) 



cos^ C, cos^ a + sin^ ji svc? a 



(2.6) 
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. sin C COS a ,„ „, 

sm{D-F) = , . . ' ^ . , (2.7) 

V sm Q COS'' a + cos^ fj. sm a 

. sin u cos Of 

sHD + F) = . / (2.8) 

V sm /i sm a + cos^ cos^ a 

The parameters, D, E and can be assumed to be in [0, 2tt) if needed. 

The next item on the list is the foUowing decomposition of any matrix, 5' 6 SU{4) into 

factors which are exponentials of /y , i = 1,2, j — x,y and the matrix Iizl2z, fl^ . 



5 = ii:i0i^2e"'''^^^''e"'^^2"e"'^i^i^-^2^e"'T'^i"e"*T'^="e"*T'^i^e"'T'^^="e~'''^'^i^-^^^e"^ 

(2.9) 



In (2.£) the matrices Ki,i = 1,...,4 are some matrices in SU{2). This decomposition 
follows from the well known fact that the Lie algebra su{2) su{2) and its orthogonal 
complement in the Lie algebra sm(4) provide a Cartan decomposition of the Lie group 
SU{A), In this fact and some calculations are used to obtain the decom- 

position ( |2.9| ). Note, however that they do not provide any formulae for the Ki,i = 1, . . . , 4 
and the Ok, k — 1, ... ,3 in terms of the target matrix S . In the appendix, we will amelio- 
rate this problem. For the moment, however, expanding each of the Ki, i = 1, . . . , 4 into 



its ax, cTy Eulcr angles via Equation (2.6) and using some Kronecker calculus, leads to the 



following equation, which is the one we will work with: 

g _ gi-Dl/lxgi-Bl/l!,gi-Fl-flxgi-D2/2xgi-B2/2j,giF2/2xg-if-ri!,g-if /2yg-i£'l-flz-f2z ^2 10) 
g-*T^-flyg-«^-f2yg-iT=^-flxg-i^/2xg-ie2-fl^-f22 

p-i^Ii:, p-i^h:^ p-iSshJ^z piDsIi^ iE^hy iFsI^^ iD2l2^ iE2l2y iF.2l2!, 



The {D„Ei,F,) are the Euler angles of the G 817(2), i = l,...,A. 



3 Determination of the Frequencies and Phases 

Consider a pair of coupled spins in the weak coupling limit. The system is assumed to be 
heteronuclear, so that it is possible to address each spin individually as long the as the 
frequency of the corresponding field is resonant with the Larmor frequency of the spin, LUi, 
in question and the field can be bounded in amplitude. Thus, the model is 
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where the internal Hamiltonian, A is wiCTz ^ /2 + '^2^2 ^ ctz + J<^z ® , and the interaction 



Hamiltonians are Bi = bidx <E) h and B2 — bidy (E) l2- The bi,i = 1,2 (62 appears below) 
are constants related to the gyromagnetic ratios, and the tOi^i — 1,2 are the Larmor 
frequencies. Finally, ui{t) and U2{t) are sinusoidal fields to be designed: 



The frequency, uj, will be taken to be toi. So the design procedure amounts to specifying 
the amplitude, c, the phase, (jj and the duration of the pulses. These will be chosen in 
a piecewise constant manner. The fourth section is essentially devoted to finding the 
amplitudes and the durations. How the phases ought to be chosen will become clear later 
in this section. 

The above equation was derived assuming that the first spin was being addressed. 
If the second spin is being addressed, then Bi,B2 would be replaced by 62/2 €5 ctx and 
&2-f2 €5 Uy respectively, while the frequency, oj, of the field would be replaced by uj2- 

Let us now derive the rotating frame in which the problem of preparation of a target 
will be translated into finding picewise constant controls for an associated system in the 
rotating frame (see |2^ for related considerations). Set 



Ui{t) — c cos{ujt + (/)), M2 = c sin(a;t + 0) 



(3.11) 



U{t) = e*^V{t) 



with 



F = 



-(^2-^2 ®(Jz+ ^iCTz (X) I2) 



Then a few calculations reveal that 



il = --{JhzhzW - -{cbiA I2)U 



(3.12) 



where 








A 



V 
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Note that the matrix A is independent of time and is parametrized by the phase, 4>, of 
the field which can he choosen. Thus, if we choose (/) = 0, A is tTa; and if we set = -I, A 

is Gy. 

A similar calculation reveals that the same rotating frame can be used to address 
the other spin (with the frequency of the field, lo — 72) to obtain the following equation 

U = -'^{JhzhzW - '-{cb2l2 ® A)U (3.13) 

Once again by choosing the phase, 4>, one can ensure that A is either a-y or cTx- 
The upshot of the foregoing is that by choosing the frequency of the field to be resonant 
with one of the spins and by choosing the phase in an appropriate manner, passage to a 
unique rotating frame leads to the following system, which is controlled by constant inputs: 

U ^ -^AU - ^dBU,U e SUiA) (3.14) 

with A = Juz ® <Jz and B one of the matrices Iix , I2X , Iiy , l2y where Iij = (Jj ® l2, j — x,y 
and l2j ^ I2 ® '^jij — ^tV- The constant, d, is related to the amplitude of the field 
and other constants of the system. This is useful because the Cartan decomposition, i.e., 



Equation (2.10) of Section 2, consists precisely of the exponentials of one of A,Iij,i — 
l,2,j — x,y. Thus, if the exponential of a certain /.y is required we choose the frequency 
and phase of the field so that B of Equation ( 4.15| ) becomes Iij, and then follow the 



procedure in the next section to determine, d (hence, the amplitude of the field). 

It is interesting to observe that while in the rotating frame, we are exciting one spin 
with one of the x or y magnetic field components, in the original frame we are exciting 
any one spin by using both the x and y magnetic field components. 

4 Determining the Amplitudes 

The basic model derived in the previous section leads to the following system in the rotating 
frame, for which we will design piecewise constant controls: 

U = -iAU - iu{t)BU, U{Q) = h (4.15) 
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where A = Iizl2z — <^z ®<^zi and B is one of the lij^i — 1,2, j — x,y and u{t) the control 
(to be determined) is piecewise constant. 

Note, for the sake of easy bookeeping, the constants J, hi and i have been dropped in 



Equation (4-15). However, once results for the above model are available it is a routine 
matter to derive results for the actual model, with these parameters present. 

Going back to the equation ( ^.10 ) it is clear that factors which are the exponentials 



of —ilizhz (i-e., of —iA) can be generated by free evolution i.e., by setting u{t) = for 
an amount of time given by the corresponding 6k- Indeed, e"*^'''^ is a diagonal matrix, 
diag(e~*^'=,e*^'=,e*^'=,e~*^'=). If 6k > 0, then e"*^*^ can be prepared by free evolution for 
9k units of time. U 6k < then, e^*^*^"^ can be prepared by free evolution for 2tt + 6 
units of time. Only for free evolution terms are periodicity arguments resorted 
to. For other factors, instead of resorting to periodicity arguments, substantial use of the 
structure of the drift ~iA will be made. 

So all that remains to be addressed is preparing factors which are exponentials of the 
—ilij via controlled pulses. In other words, the main goal at this stage is to decompose 
the exponential, e~*^^'^ , i = l,2,j — x,y, for any L <E R, as: 



-iLIi- 

e 



n^^ie(-^'^'=^-^'"=^-),z = l,2,j=x,y, A = hJ2z (4.16) 



satsifying i) Ol ak > 0, fc = 1, . . . ,Q and ii) 02 | bk |< C,k — 1, . . . ,Q. Constructive 
methods for obtaining the desired ak, bk in terms of S will be provided. At this stage the 
formulae for the Ok and bk will be provided in terms of the parameters, (Di, Ei, Fi),i — 



1 , . . . , 4 and 6k, k — I, ... ,3 of Equation ( 2 . 1 C ) The appendix indicates a procedure which 



obtains these S'C/(4) "Euler angles" directly from the entries of S. Note that ak is the 
duration of the kth pulse and ^ is its amplitude. The basic strategy is to first achieve 
Ol and then constructively modify the resulting decomposition to meet 02 also. The 
proof below shows how achieving 02 leads also to arbitrarily bounding the amplitude of 
the pulse. 

Theorem 4.1 The matrix exp{—iLIix), L ^ R can be factored explicitly as H'l^^cxp{—iakA' 
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ibkhx) with A = hzhz and with afe > 0, fc = 1, . . . , 3 and 6i = = 63, if cosL 7^ 0. If 
cosL = 0, then four factors are needed, with bi = = 64. Further, this factoriza- 
tion can be refined constructively, by increasing the number of factors, to ensure that 
\ bk \< C,k = 1, . . . ,Q. This decomposition can be further refined to meet the condition, 
I ^ |< D, for any prescribed D > 0. Similar statements hold for the exponentials of 
hxjhy and l2y with the matrix for appropriately different sets of a^'s and bk- 

Proof: The proof will be given only for the exponential of exp{—iLIix), since the 
proof for the others are very similar. 

The special structure of hzhz and Ii^ results in the following matrix for exp{—iakA — 

ibkhx) 



I 



cosAfe - j^ak sinAfe 



■ sin At 



V 



COS A/c + j^ak sin Afc 




sin Afe 



-^SinAfe 



cosAfe + f-Ofc sinAfe 







-ibk 



sin Afe 



cos ^k — j^dk sin Afe 
(4.17) 



where Afe = ^ (a'f. + b'f,. In particular, setting afe = and bk = L yields the exponential 
of —iLIix- Similarly, an explicit calculation yields the following formula: 



exp{—iL(7y (g) az) 



V 



cosi 





sinL 








cosL 





— sinZ/ 


— sinL 





cosL 








sini 





cosi 



Using the last two formulae, the following useful identity is obtained: 

exp{-iLIix) = exp{-i—Iizl2z)exp{-iLay (g) cTz)exp{-i-Iizl2z) 



(4.18) 



(4.19) 



Thus, to prepare the matrix exp{—iLIix), free evolution for ^ units of time for the first 
factor above is used and free evolution for j units of time for the third factor is used. 
Therefore, it remains to produce the middle factor, exp{—iL<jy®(Jz) via controlled pulses. 
Since, the matrix —iLoy ® Oz is, upto a constant, the commutator of A and -Bi, it seems 
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plausible that its exponential can be represented as a product, Il^^iexp[—i{aklizl2z + 

bkhx)]- 

To demonstrate this, three cases need to be considered: i) cosL > 0, ii) cosL < and 
iii) cosi = 0. 

The Case cosL > 0: Evaluating the product Il1^-^cxp[—i{aklizl2z + bklix)] and choosing 
Ai = \Ja\ + b\ = ^ and A2 = ■\/af+6| = J and equating the result to the matrix 
exp{—iL<jy <Jz) leads to the following equations: 

a\a2 + 6162 

—7 = cos L 

A1A2 

ai&2 - a,2bi . r 

= -""^ 

These two equations can be solved as follows. Choose 61 = 0, oi = 62 = — f sinL, 02 = 
f cosL. Thus, 01 has been met for the case that cos L > 0. 

The Ccise cosi < 0: Now choose Ai = | and A2 = f . The resulting set of equations has 
the following solution: 61 = — f sin L,ai = — ^cos i,a2 = f and 62 = 0. Thus, Ol has 
been achieved for this case also. 

The Case cosi = 0: Now choose Ai = | and A2 = f . Then, the equations to solve 
become °^'\^+^;''^ = and = (-1) sin(i). If sin(L) = 1, choose ai = ^Ai and 

a2 = -^^2, and 61 = — ai and 62 = 02- If sinL = —1, then choose ai = "^Ai, a2 = "^-^2, 
and 61 = ai and 62 = — ^2, 

Now concatenating the pulses which prepare exp{—iLay (g) az) with the free evolution 
terms which prepare exp{—i^az cTz) and exp{—ijaz <7z) yields the stated values for 
the number of factors for Ol. 

To meet 02, notice only the e"'^'^"®'^^ term needs to be addressed, since the others 
are free evolution terms and hence have pulse area equal to 0. Even in the preparation 
of the Q-i'^'^v^^" term, there is only one term which is not a free evolution term (except 
when cosi = 0). The corresponding, | 6fc | is exactly equal to ^ | sinL |. So to meet 
02, 6 has to be such that | sinL |< C^. This amounts to requiring that L be within a 
prescribed bound of 0. If L is not already of the form, then we factor e'^^^^" as 11^^ ■^e*^''^^^ 
with the Lfe satisfying the required deviation from condition. Clearly this can always be 
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done. 

Notice further, that this process of meeting 02 also ensures that the amphtude of the 
pulse, I ^ I can also be bounded arbitrarily. Indeed, the amplitude of the pulses are either 
(corresponding to free evolution terms) or | tanL |. Clearly any process which ensures 
that I sin L \ is within a prescribed bound can ensure the same for | tan L \ . 

If cos L — 0, then write exp{—iLIix) as {exp{~i^Iix))'^ , and proceed as in the cosL 7^ 
cases. This finishes the construction. 

Remark 4.1 The values obtained for at and bk are certainly not the only possibilities. 
Since one of the principal goals of this paper is to show that the pulse amplitudes can 
be arbitrarily bounded, the and bk satsifying Ol were so chosen that the resulting 
decomposition could be modified with minimal fuss, to meet 02. This means that the 
proof chosen was biased towards free evolution terms. In practice, of course one can find 
other values so that the cumulative time taken can also be kept within reasonable bounds. 

Remark 4.2 Suppose S is prepared in the rotating frame in T5 units of time. Then in 
the orginal coordinates the matrix e^'^^^ S has been prepared. Depending on the system, 
it may be desirable to rectify this error. Since, Ts > [with 6k given in Equation 

( |1.2| )], this problem cannot be wished away by hard pulses. Of course, hard pulses can be 
further used to generate e^^^ . However, this introduces further inaccuracies in addition to 
those caused by the use of hard pulses to prepare S in the rotating frame. We suggest two 
methods which do not need hard pulses to rectify this deviation. The first is to prepare 
gToF g gQj:^ pulses in the rotating frame for a real parameter. To, to be chosen such that 
the time, Ti, to prepare e^"^ in the rotating frame, satisfies Ti +Ts = To. This leads to 
a transcendental equation for Tq [ specifically, + -^)^ + Ts] = To — ^(coswiro + 

cos(jj2To)]. This has to be solved numerically. The second method is to use optimal control 
to drive the system, in the original coordinates, from I4 to e^^^ , with a quadratic cost 
functional incorporating bounds on the field fluence and deviation of the state from the 
target. The fact that the target state, e^-^''^ is diagonal will help reduce the complexity of 
optimal control calculations. 
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5 Conclusions 

In this paper, a constructive procedure for generating a desired unitary generator in SU (4), 
via the control of a coupled, heteronuclear, two spin system was described. The sinusoidal 
pulses that were produced were not hard pulses, but instead could be bounded both in 
amplitude and pulse area. 

An interesting problem is to investigate the methodology of this paper for systems 
which are not studied by addressing individual spins selectively. It is relatively straight- 
forward to see that an analogous Cartan decomposition can be modified to express every 
target 5 in a manner analogous to Equation ( |1.2| ). However, the remaining calculations 
seem to require new methods. Investigation of this problem will be a worthwhile research 
problem. 



6 Appendix: Determining the Euler Parameters for 

SU{A) 

To make the methodology proposed here (or for that matter any methodology based on the 
Cartan decompositions) genuinely constructive, it is extremely desirable to find the fifteen 



parameters, {Di, Ei, Fi),i = 1, . . . , 4 and 9k, k — 1, . . . , 3 of Equation (2.10), as explicitly as 
possible, in terms of the entries of the target matrix S G SU{A). In principle, this amounts 
to solving sixteen equations for fifteen unknowns. However, this is not really a satsifactory 
state of affairs. To make the point further clear, consider determining the Euler angles of 
an SU (2) matrix in closed form in terms of the entries of the matrix. This is principally 



facilitated by the availability of the Cayley-Klein representation, i.e.. Equation (2.4). A 
similar representation for 5t/(4) matrices is not available in the literature, and thus the 
resulting system of equations cannot even be written down in a manner which will facilitate 
investigating the possibility of a closed form determination of the Di, Ei, Fi,i — 1, . . . , 4 
and the 6k, k = 1, . . . , 3. 

Therefore, we eschew working directly with S itself. Instead, S will be first decomposed 



14 



into a product of matrices, each of which has a simpler structure, so that finding these 
fifteen parameters for each of these factors is more tractable. Note that strictly speaking. 



this will not result in a factorization of the form in Equation (2.10) for the original matrix 
S, since now the factors of each of the Si become intertwined. However, it will produce a 
factorization of the form in Equation (pT^). But, this is all that is needed to achieve the 
constructive generation of S. 

The main idea is to write S = H^-^S', so that each of the Si can be expressed solely 
in terms of a single SU{2) matrix. While, this may mean that the number of parameters 
(and, thus the number of control pulses) needed are larger than what would result if the 
fifteen parameters were directly determined for S, it has the advantage that close to explicit 
formulae can be produced for the Euler angles, whereas (pending further investigation) 
there is nothing remotely close to explicit when working with S itself. Furthermore, as 
will be clear soon, many of the Kk, k — 1, . . . , 4 turn out to be I2 for many of the factors 
5i, and likewise many of the 9k are zero. Thus, the number of factors and pulses is not all 
that high. 

The decomposition that will be used is the standard Givens decomposition, 
modified slightly for the problem at hand. The usual Givens decomposition is produced as 
follows. Premultiply S\ the inverse of 5, by a sequuence of matrices, S'i, which successively 
reduce the columns of S to the unit vectors, ei, i.e., the fourth column is reduced to 
(0, 0, 0, 1) and the third to (0, 0, 1, 0) etc.. Then it can be shown that Q — [||, and thus 

s = nj_2 Si 

Usually these Si are taken to be a matrix, which upto permutation of rows and columns, is 
a block matrix consisting of I2 and an explicitly determined SU{2) matrix. Here, we will 
make a slight modification, we will take 5*6 to be a tensor product of two SU{2) matrices. 
The remaining five will be, upto permutation, block matrices with blocks equal to I2 and 
a specific SU{2) matrix. 

To briefiy illustrate the structure of 5*6 (the remaining Sk are constructed in the manner 
descirbed in [^), suppose that Suppose S^ — col{a,b,c,d) (thus a,b,c,d are the four 
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columns of S^). 

Then choose 5*6 = e'^'^" (g) S{aQ, Ce, Me) where S{aQ, Ce, Me) is the unique SU (2) matrix 
which takes the vector (^1,^2) to the vector (|| (^1,^2) ||,0). Note that once di,d2 are 
known one can exphcitly write down this SU{2) matrix (here use is being made of the fact 
that given any two points on a sphere of any radius in there is a unique SU (2) matrix 
which conveys the first to the second and that this matrix can be found explicitly). 



The remaining, Sk have the following structure 



and 



1 
1 

cos as e'^^ sin 0:56 













V5 



sinase'^'^ cosase ^ 



cos 0:46'^* 



sin 0:46 

1 
1 



sinQ!4e'(''-'**) cos 046"'^* ^ 



cos as e'^^ sin as e 

sinaae*^''-''^) cosase-*'^^ q 






e'l"^ 



1 
1 















COS a2e 



<2 



sin 026*''^ 



sma2e'(''-''2) cosa2e-''^2 



V 
/ 



51 = 







cosaie'fi sinaie''^! 
cosaie-'«i 
1 







sinaie 













1 
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Once again each of the SU{2) matrices, S'(afc, Cfe, ^fe), k = 1, . . . , 5, can be construc- 
tively determined from the entries of S. 

Now let us determine the fifteen Euler angles for the matrices 5'i, i = 1, . . . , 6. Since, 
the Si are significantly simpler than the matrix, 5, this is a more tractable task. In 
what follows the matrices Kk,i = 1, . . . , 4 and the real constants, 9k, k — 1, ... ,3 will be 
described for each of the Si,i = 1, . . . , 6. Once the are known, it is easy to find the 



fJxT^y Euler angles, {D^, E^, F^), by using Equation (2.6). Therefore, that step will not 
be executed here. 

To simplify notation, three real parameters P, Q and R are introduced. They are 
related to the parameters 0^, A; = 1, . . . , 3 by: 

P^^l_^, Q^^, R^h±h (6.20) 
4 4 4 

Below, the values of P, Q and R for the Si,i — 1, . . . , 6 will be given. Obtaining the 

9k, k = 1, . . . , 3 is then routine. 

Sq: Clearly, Sq = Ki®K2, where Ki = e'^^'y and K2 = e'-°eT.g-iE6a„g-iF6<T,^ rpj^g g^jg^. 
angles {Dq, Eq, Fq) are the Euler angles of the matrix S{ae, Ce, fJ'e) determined according 
to Equation (^^. The matrices K3, K4 can be taken to be I2 and the constants, P, Q and 
R (and hence the 9k, k — 1,2, 3) can be set equal to zero. 

SqjSi and 53: The matrices Si and S3 are essentially the same in structure, and they 
are analogous to S5. Thus, calculations for ^5 will be shown here and the modifications 
required for Si and ^3 will be given. 

So consider determing the fifteen Euler parameters of 6*5. 

Pick Ki = I2 and K2 ~ S{a,(^,fi) where the parameters (a,C,/i) will be presently 
determined. Choose P an i? equal to zero (thus ai = = 02). Multiplying out all but 
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the factors ® K4, leads to the foUowing matrix: 



cosae*('^-'3) sinae*('^+'3) 



cos ae 



i(Q+0 



sm ae 






»(m-Q) 



V 



sinQ;e^(''-''+'3) cos ae-^^^+O 



sinae*(''"^"'5) cos ae-'^'^-Cs) 


^ 

Now choose = I2 and K4 to be the inverse of the top left hand block of the last matrix 
(note since the matrix in question is in SU (2) it is very straightforward to find its inverse). 

This then means that the matrix ^(as, Cs, Ms) (which is known) should equal the ma- 
trix: 



2 ^^i2Q 



COS ae 



+ sin'^e''^^ cos a sin a(e*(2Q+A'+C-^) + g-*(M+C-2Q) ) 



sin^ ae'^^^ + cos^ ae 



^ cosasina(e*(''-'"-«+2Q) _,_g-*2Q 

We will find Q!,CjM ^^'^ Q (and hence as = AQ) by equating the top row of the two 
matrices in the last equation (since both the matrices are in SU{2), this will automatically 
mean that the second row of the two are the same). This leads to the equations: 



cos^ 2Q + cos^ a sin^ 2Q — cos as 
cos 2a tan 2Q — tan ^5 
tan /j,5 



^ cos(m + C) ^ 



sin(^ + C) ■ 

In the above system of equations, the unknowns are a, C, fi and Q, and the known variables 
are asiCsiMs- Thus, we have three transcendental equations for four unknowns. Thus, 
there will be in general many solutions. Notice, that the first two equations involve only 
two of the unknowns, viz., a and Q, and thus it is this pair of equations which will have 
to be solved numerically. The one parameter family of freedom comes from the third 
equation, where there are two unknowns. 

In summary, the Euler angles of S5 have been determined. The Euler angles of S3 
(and, thus Si) can also be determined via a similar technique. Indeed, the calculations are 
similar to those for 5*3 except that K4 will be taken to be the inverse of the SU{2) matrix 
in the bottom block of the matrix resulting from multiplying all but the last six factors of 



Equation (2.1C) 
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S2 and 5*4: The matrices S'4 and S2 bear a resemblance to one another. Calculations for 
S2 will be shown together with the modifications needed for the S'4 case. 

Choose Ki = e"'!'"^, K2 — e"'^'^^, K3 = I2 and K4 = e*''^'^^, for some real numbers 
rik,k — 1, . . . , 3 to be determined shortly. Choose Q = and P — and R — a2 (where 
a2 is a Cayley-Klein parameter of 5*2 = S{a2, C2, ^2), and thus, is known beforehand). 

The parameters 77^, A: = 1, . . . , 3 are found by solving the linear system of equations: 

m = 
V1-V2 + m = C2 

TT 

m-m-m = m + ^ 

For S4 the only modification needed in this procedure is that R = and P = a^. 

This completes the determination of the fifteen, SU{4:), "Euler" angles for each of the 
matrices Si,i = 1, . . . , 6 in the specially chosen Givens decomposition of the given target 
S e SU{4). 

Remark 6.1 There is considerable liberty in the Givens decomposition. For instance, 
the order in which the columns are reduced to the corresponding unit vectors is one such 
degree of freedom. The factors that were chosen above were expressly intended to facilitate 
the calculation of the corresponding fifteen St/ (4) angles. It is an interesting problem to 



find other factorizations which yield different values of the tk in Equation (1.2) 
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